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Quantum systems with discrete symmetries can usually be desymmetrized, but this strategy 
fails when considering transport in open systems with a symmetry that maps different openings 
onto each other. We investigate the joint probability density of transmission eigenvalues for such 
systems in random-matrix theory. In the orthogonal symmetry class we show that the eigenvalue 
statistics manifests level repulsion between only every second transmission eigenvalue. This finds its 
natural statistical interpretation as a staggered superposition of two eigenvalue sequences. For a large 
number of channels, the statistics for a system with a lead-transposing symmetry approaches that of 
a superposition of two uncorrelated sets of eigenvalues as in systems with a lead-preserving symmetry 
(which can be desymmetrized). These predictions are confirmed by numerical computations of the 
transmission-eigenvalue spacing distribution for quantum billiards and for the open kicked rotator. 

PACS numbers: 05.45.Mt, 05.60.Gg, 73.23.-b 



I. INTRODUCTION 

Mesoscopic systems exhibit variations in their phase- 
coherent electronic transport properties that are con- 
veniently characterized via statistical approaches. Ge- 
ometries that classically give rise to chaotic motion typ- 
ically display universal fluctuations which can be cap- 
tured using ensembles of random scattering matrices llj . 
For normal conductors the universal properties fall into 
Dyson's three universality classes with symmetry index 
/3 = 1,2,4 0, while a further seven universality classes 
can be identified in the presence of superconducting or 
chiral particle- hole symmetries [1]. A powerful tool to 
distinguish these ensembles is the amount of level re- 
pulsion between the transmission eigenvalues Tn- These 
eigenvalues determine fundamental transport properties 
such as the conductance G or the shot-noise Fano fac- 
tor P, In the Dyson ensembles, the probability 
density to find two closely spaced adjacent transmission 
eigenvalues with small distance s = Tn+i — T„ is sup- 
pressed as P{s) oc s'^ 5, 6]. This introduces a stiffness 
in the transmission-eigenvalue sequence which suppresses 
the fluctuations of the conductance and of the Fano factor 
when compared to the case of uncorrelated transmission 
eigenvalues (the latter being characteristic for classically 
integrable systems with a complete set of good quantum 
numbers) [ll, 

From the investigation of closed systems it is well 
known that discrete symmetries result in a reduction of 
level repulsion. In such systems, desymmetrization deliv- 
ers independent variants of the system which differ by the 
boundary conditions on the symmetry lines (e.g., Dirich- 
let and Neumann boundary conditions for eigenfunctions 
of odd and even parity, respectively). The statistics of the 
desymmetrized versions can depend on the dimensional- 
ity of the irreducible representation Q , but still remain 



within the conventional universality classes. The com- 
bined level statistics is then built by superimposing the 
independent level sequences of the desymmetrized vari- 
ants @. In open systems, this concept of desymmetriza- 
tion can be directly applied as long as the symmetry in 
question preserves the shape and position of the leads 
li. 

This paper is motivated by the observation that sys- 
tems with a lead-transposing symmetry (which maps dif- 
ferent openings onto each other while leaving the dynam- 
ics in the system unchanged) exhibit transport properties 
that can only be understood as collective features of the 
desymmetrized variants of the system Q. An obvi- 
ous indication of this complication is the fact that the 
symmetry-reduced variants only possess a reduced num- 
ber of leads (we concentrate on systems with two leads, 
for which the desymmetrized variants only possess a sin- 
gle lead). We demonstrate that such systems exhibit nev- 
ertheless a reduced repulsion of transmission eigenvalues 
which is similar to that for systems with a lead-preserving 
symmetry. For a large number N of transport chan- 
nels, the local statistical fluctuations in the eigenvalue 
sequence indeed become indistinguishable for both types 
of symmetry. However, for a small numbers of channels, 
the statistics differ from each other, which can be traced 
back to the absence or presence of 1/iV corrections in 
these ensembles. 

In the specific case of /3 = 1, we derive exact closed ex- 
pressions for the joint probability density of transmission 
eigenvalues thereby gaining detailed insight into these 
statistical features. In particular, we find for both the 
lead-preserving and the lead-transposing symmetry class 
that level repulsion occurs only between every second 
transmission eigenvalue. The fluctuations in the trans- 
mission eigenvalue sequence hence find their most natu- 
ral statistical interpretation in a staggered superposition 



of two independent level sequences. In such a superposi- 
tion, the transmission eigenvalues alternate between the 
two sequences when they are ordered by magnitude. 

The exact expressions for the joint probability density 
with (3 = 1 are different for the two types of symme- 
try. Hence, the details of the transport statistics for a 
lead-transposing symmetry deviate from those for a lead- 
preserving symmetry. We show that these deviations are 
most significant for a small number of channels, while for 
a large number of channels the local eigenvalue statistics 
do indeed converge onto each other. 

Previous studies of open systems with lead-transposing 
or lead-preserving symmetries have derived the distribu- 
tion of transmission eigenvalues for one or two open chan- 
nels and the one-point density for arbitrary numbers of 
channels [1, [^, [13, [lH ■ For time- reversal symmetric sys- 
tems with P = 1, a key observation of these works was an 
enhancement of universal fluctuations for both types of 
symmetry (when compared to asymmetric systems). For 
systems with a lead-transposing symmetry it was found 
that the weak localization correction is vanishing, lead- 
ing to ensemble averaged expressions for the conductance 
and for the shot noise Fano factor which are entirely in- 
dependent of the channel number N . The underlying 
staggered level statistics embodied in the joint distribu- 
tion of transmission eigenvalues provides a unifying ex- 
planation for all of these observations. We verify our 
predictions by numerical computations for quantum bil- 
liards [H, [3 and for the open kicked rotator Ti', IT^, [l^ . 

This paper is organized as follows. Section [Til provides 
background information on the scattering approach to 
transport and on standard random-matrix theory. In Sec. 
mil we revisit the case of systems with a lead-preserving 
symmetry and provide the exact reformulation of the 
eigenvalue statistics in the orthogonal symmetry class 
(/3 = 1) as a staggered superposition of two eigenvalue 
sequences. Section IIVI concerns systems with a lead- 
transposing symmetry. In particular, for (3=1 we derive 
the exact joint probability density of transmission eigen- 
values for arbitrary iV, and show that this again takes the 
form of a staggered eigenvalue sequence. We also describe 
the convergence of the local statistics for both types of 
symmetry, which emerges in the limit N —^ oo. Section 
fVl provides numerical results that illustrate the similari- 
ties and differences of the random-matrix ensembles for 
the two symmetry classes. This section also contains the 
comparison to specific model systems. Section I VII pro- 
vides a summary and discussion of our main results. 



II. BASIC CONCEPTS 

A. Scattering approach to transport 

Figure [T] depicts open two-dimensional quantum bil- 
liards representing mesoscopic systems with two attached 
leads (L - left and R - right), each carrying N in- 
coming and N outgoing modes. The systems in Fig. 





(d) 




FIG. 1: (Color online) Sketches of quantum billiards (a,b) 
without any spatial symmetry, (c) with a lead-preserving 
symmetry, (d) with both a lead-preserving as well as a lead- 
transposing symmetry, (e) with a lead-transposing reflection 
symmetry, and (f) with a lead-transposing inversion symme- 
try. The inversion symmetry in panel (f) survives in the pres- 
ence of a finite magnetic field, as is indicated by a symmetric 
pair of trajectories. 



[Ha,b) are asymmetric while those in Fig. me,f) possess 
a lead-transposing reflection or inversion symmetry, re- 
spectively. The inversion symmetry in panel (f) survives 
in the presence of a finite magnetic field, which how- 
ever breaks time-reversal symmetry (an inversion symme- 
try does not induce a generalized antiunitary symmetry, 
in contrast to refiection symmetries [gI]). In the middle 
panels. Fig. [Ijc) shows a system with a lead-preserving 
symmetry, and Fig. [Ijd) shows a system which possesses 
both a lead-preserving (up-down) and a lead-transposing 
(right-left) symmetry. 

In order to describe the phase-coherent transport 
through these systems for small bias voltage V, one solves 
the Schrodinger equation for fixed values of the 27V am- 
plitudes a = [ai^\ a^^]"^ in the incoming modes. This 
results in linear relations b = Sa for the 27V amplitudes 



27V X 27V-dimensional scattering matrix of the form 



b = [bn', blC'Y the outgoing modes, which delivers a 



5* 



r t' 
t r' 



(1) 



Here r, r', i, t' are 7V x 7V-dimensional matrices describing 
reflection at each lead and transmission from one lead to 
the other, respectively. 

The scattering matrix is unitary, and its structure is 
further constrained by symmetries of the system. The 
three main universality classes arise for systems with 
time-reversal and spin-rotation symmetry (orthogonal 
symmetry class with S = S'^ , symmetry index /3 = 1), 
systems without time-reversal symmetry (unitary sym- 
metry class with no constraints on S, (3 = 2), and sys- 



3 



terns with time-reversal but broken spin-rotation sym- 
metry (symplectic symmetry class composed of self-dual 
matrices S = S^, (3 — 4). Spatial symmetries entail ad- 
ditional constraints on the scattering matrix, which are 
detailed in Sees. |TTT] and |TVl 

The transmission eigenvalues T!„ are defined as the 
eigenvalues of the hermitian matrix tt'' . In the case of 
spin-independent transport or Kramers degeneracy (the 
latter occurs for P — A), the transmission eigenvalues are 
twofold degenerate. We then only account for each pair 
of eigenvalues once and introduce a spin-degeneracy fac- 
tor a = 2. When the two- fold degeneracy is lifted then 
a = 1. From here on, N refers to the number of distinct 
transmission eigenvalues (ignoring accidental degenera- 
cies). Furthermore we will assume that the transmission 
eigenvalues are ordered by magnitude. 



Ti < To < T, < 



<Tn, 



(2) 



as this results in a number of technical simplifications. 
The conductance quantum is defined as Go = ae^ /h. 

With these conventions, the transmission eigenvalues 
determine fundamental transport properties such as the 
conductance via 



G — Go T„ 



(3) 



and ultimately results in conductance fluctuations of the 
order of a single conductance quantum, which for large 
N approach the asymptotic value 



varG/Go 



1 

8^' 



(6) 



For large N, the one-point probability density of trans- 
mission eigenvalues approaches 



P{T) 



1 



(7) 



The second product in Eq. ([5]) induces an asymmetry into 
this bi-modal distribution, which for large TV results in 
the weak-localization correction 



(G) - |go - Go (J - ^ 



(8) 



of the ensemble-averaged conductance. 

An insightful quantity derived from the joint proba- 
bility density P({r„}) is the distribution P{s) of spac- 
ings s = Tn+i — Tn between neighboring transmission 
eigenvalues. For uncorrelated eigenvalues with average 
spacing s one would expect a Poisson distribution. 



and the shot-noise power via 



N 



P = 2Goey^r„(l-T„). 



(4) 



n=l 



Here V is the bias voltage, which is assumed to be small. 



B. Dyson's circular ensembles 

Random-matrix theory delivers a statistical descrip- 
tion of transport by drawing the scattering matrices from 
ensembles of unitary matrices which obey the constraints 
of the given universality class. For the three main univer- 
sality classes with /? = 1, 2, or 4, random-matrix theory 
is based on Dyson's circular ensembles, for which the 
probability measure is given by the Haar measure of uni- 
tary symmetric, unitary, or unitary self-dual matrices, 
respectively. The joint probability density of transmis- 
sion eigenvalues then takes the form p| 



(5) 



The first product in Eq. ([5]) involves pairs of transmis- 
sion eigenvalues and favors sequences in which neighbor- 
ing transmission eigenvalues do not approach each other 
closely. (As we have ordered the transmission eigenval- 
ues by magnitude, all differences Tm — Tn are positive.) 
This suppresses fluctuations in the eigenvalue sequence 



Pis) 



(9) 



while for the circular ensembles and A'' ^ 1, the spac- 
ing distribution can be well approximated by the Wigner 
distributions d, @, 



Pis) 



2^sexpi-||3 



32 



exp 



36^5 exp 



i-m 

I 64s' 



(3 = 1 
(3 = 2 

(3 = 4: 



(10) 



Lead-preserving and lead-transposing symmetries en- 
tail further constraints on the scattering matrix. The 
consequences of these constraints for the transmission 
eigenvalue statistics are explored in the remainder of this 
paper. 



III. LEAD-PRESERVING SYMMETRIES 

A useful reference point for our subsequent investiga- 
tion of systems with a lead-transposing symmetry (in Sec. 
IIV[) are open systems with a lead-preserving symmetry, 
to which one can directly apply the standard ideas of 
desymmetrization. The goal of the present section is to 
reformulate the resulting random-matrix statistics for the 
case of a lead-preserving symmetry with /3 = 1 as a stag- 
gered level repulsion, as this will allow us to establish a 
connection to the case of a lead-transposing symmetry. 
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B. Conventional random- matrix theory 
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FIG. 2: (Color online) (a,b) Desymmetrization of the quan- 
tum billiard with a lead-preserving reflection symmetry, 
shown in Fig. [TJc). (c) Sketch of the individual transmis- 
sion eigenvalue sequences of fixed parity, (d) Reorganization 
as a staggered level sequence, where transmission eigenvalues 
alternate after ordering by magnitude. 



Within random-matrix theory, the joint probabihty 
distribution of the total transmission-eigenvalue sequence 
is the sum of the corresponding probabilities for each way 
to distribute the transmission eigenvalues into two sets 
containing A^i and N2 eigenvalues. With each sequence 
obeying the statistics of the appropriate Dyson ensemble 
one finds with Eq. ^ 

F({r„}) cx ^ n iTi^-T.j' n {Tj^-Tjy 

>n ■m'>n 

\{t-^^^'\ (11) 



/eP m>n 

N 

X 

1=1 



A. Constraints on the scattering matrix 



For large N , the separation into two effectively inde- 
pendent systems with Dirichlet and Neumann boundary 
conditions naturally results in a doubling of the con- 
ductance fluctuations ([6]) and a doubling of the weak- 
localization correction ([5]). Moreover, level repulsion is 
only effective for transmission eigenvalues which are part 
of the same sequence. This modifies the spacing proba- 
bility density, which can be calculated from the general 
expression [Sl] 



An example of a system with a lead-preserving reflec- 
tion symmetry is shown in Fig. [1] (c). Figure [5] (a,b) 
shows the desymmetrized version of the system, which 
is halved at the symmetry line. Dirichlet boundary con- 
ditions on the line of symmetry select scattering wave 
functions with an odd parity, while Neumann boundary 
conditions yield even parity. Consequently, the transmis- 
sion matrix t assumes a block structure where each block 
corresponds to a given parity. As dictated by the one- 
dimensional transverse-mode quantization in the leads, 
the block of even parity has dimension A^i = [(A^+ l)/2], 
while the block of odd parity has dimension N2 = [N/2] 
(here [•] denotes the integer part of a number). Hence, 
both blocks have either the same size (when N ~ N1+N2 
is even), or the block with even parity is by one larger 
than the block with odd parity (when N is odd). 

The total transmission-eigenvalue sequence is therefore 
obtained from a superposition of two sequences of size 
A'^i and N2 [for illustration see Fig. [UJc)]. In order to 
fix the way we address the elements of this superposi- 
tion, we impose the ordering of Eq. ^ and denote by 
V the set of all strictly increasing sequences of indices 
/„ e {1,2,3,..., N}, where each sequence is of length 
A^i . Such sequences are of the form I = (Ii, I2, ■ ■ ■ , Ini ) , 
where 1 < /i < /2 < /a < . . . < /jvi < N. For 
each sequence we also define a complementary sequence 
/ = (Ji, /2, Is, . . . , /atj), which consists of the indices 
1 < /i < J2 < /a < . . . < < N not contained in 
/. This partition delivers two ordered subsequences T/^ 
and Tj . 



j2 roo foo 



s + y + z)dydz (12) 



for multiple sequences i, where pi (s) is the spacing proba- 
bility densities of each sequence, while ^ is the associated 
fractional eigenvalue density. 

For two sequences following the Wigner distribution 
pO| . the resulting spacing probability densities is (s = 1) 



^xe-^'8{xl x = ^(13a) 



2 



2:i^e-^£:(x) + 



2 ' 
(13b) 



30F 



+ ^(9 + 28a;2-f 8x4)e-2-^ 
97r 



4s 

X — - — (13c) 



3V^' 



where £{x) — erfc (a;) denotes the complementary error 
function. 



C. /3 = 1: Reformulation as a staggered eigenvalue 
sequence 

In most situations encountered in random-matrix the- 
ory, the combinatorial sum over partitions involved in the 
superposition of eigenvalue sequences ([as in Eq. (jlip ] 
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cannot be performed explicitly. For the specific case 
/3 = 1, however, the combinatorial sum over / in Eq. 
(fTTj) can be carried out (see below), which then yields a 
closed-form expression 



1 



both odd both even 

. . (14) 

(A similar simplification does not present itself in the 
cases P = 2 and (3 = 4.) This result finds its nat- 
ural statistical interpretation as a staggered superposi- 
tion of two sequences, which is illustrated in Fig. ^d). 
In such a superposition, the transmission eigenvalues in 
each sequence are not distinguished by the parity of 
the associated wavefunction under the symmetry oper- 
ation. Instead, the transmission eigenvalues are ordered 
by magnitude (irrespective of parity), and one sequence 
is composed of all odd-indexed transmission eigenvalues 
(of which there are A^i) while the other sequence is com- 
posed of all even-indexed transmission eigenvalues (of 
which there are N2). Compared to the original super- 
position of two independent sequences, this differs by the 
additional constraint 



(15) 



(which is satisfied when all the ordered indices /„ are odd 
while the indices /„ are all even). 

In order to demonstrate the equivalence of Eq. (fTTj) 
(for (3 = 1) and Eq. we have to show that the level- 
repulsion terms are proportional to each other (both ex- 
pressions share the same product of one-point weights 



TliTi 



-1/2 



and the proportionality constant is fixed by 



normalization). We set out to work towards this goal by 
defining a matrix 



M 



-Vi V2 -V3 V4 -V5 Vg 
Wl W2 W3 W4 W5 Wg 



(-l)^v^ 



Wat 



which is composed of column vectors 



n T t2 rpNi~l\T 



w„ = (i,r„,r, 



rpN2- 



1\T 



(16) 

(17) 
(18) 



The determinant det M can be evaluated in two dif- 
ferent ways. In the first way, we expand it in terms of 
subdeterminants with A'^i vectors v„ from the first A^i 
rows and N2 vectors w,„ from the remaining rows. In 
other words, we sum over all determinants of the form 



neighboring rows to bring all vectors v/^ to the left 
(into row n). This results in an additional sign factor 
(_l)ai-i)+(/2-2)+...+(7„i-iVi) rpj^g determinant of the 

resulting block matrix factorizes. Overall, this expansion 
yields 



detM = (-1) 



N^{N^ + l)/2 



^ det(v/, , V/, , . . . , v/„^ ) det(wj^ , wj^ , . . . , wj^J. 



lev 



(20) 



Each subdeterminant is of the form of a Vandermonde 
determinant, and therefore 



detM - (_l)A^i(A^i+l)/2 

X E n (^^.^ 

lev m>n 



(21) 



Secondly, the determinant det M can be evaluated by 
adding in Eq. the first N2 rows to the last N2 rows. 
This yields 



det M = det 



-Vi V2 

2w2 



-V3 V4 
2^Ki 



-V5 Vg 

2wg 



Proceeding again with the evaluation of subdeterminants 
we are left with a single choice, namely, to select vectors 
v„ with odd index and vectors w„ with even index. Ac- 
counting for all signs and now also factors of two, this 
results in 

detM - (_i)^iW+i)/2 2^^ 

X det(vi, V3, V5, . . .) det(vif2, W4, Wg, . . .). (23) 

As this again involves Vandermonde determinants, we 
find 

detAf = (_i)^iW+i)/2 2^^ 

X n (Trn-Tr,) J] - Xa) ■ (24) 



The two results Eqs. (PT|) and 
able identity 



deliver the remark- 



lev m>n 



= 2 



^2 



n (Trn-Tn) J] (Trn - T„) . (25) 



det 



-Vi V2 V4 

W3 W5 Wg 



(19) 



etc., where the indices of the vectors v/^ form an or- 
dered subsequence / and the indices of the vectors wj^ 
are given by the complementary subsequence /. The al- 
ternating signs in front of the vectors vj^ can be pulled 
out of the determinant at the cost of an overall fac- 
tor (— 1)^i+^2-i-- --I-/ni ^ Next, we use permutations of 



An equivalent identity has been derived for superposi- 
tions of energy eigenvalue sequences with a length differ- 
ence of at most one which are distributed according to the 
Gaussian orthogonal ensemble 17]. Relation ((25l) shows 
that the level-repulsion term in Eq. (Ilip is indeed propor- 
tional to the level- repulsion term in Eq. . As already 

— 1/2 

mentioned, the one-point product J^^ in both ex- 

pressions is identical, and the proportionality constant is 
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(b) 



(c) 



(d) 
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FIG. 3: (Color online) (a,b) Desymmetrization of the quan- 
tum billiard with a lead-transposing reflection symmetry, 
shown in Fig. [He). (c,d) Eigenphases On of the matrix Q 
on the unit circle, and their projection Eq. ((28} which deliv- 
ers the transmission eigenvalues. 



fixed by normalization. It follows that for /3 = 1, the in- 
dependent superposition of two transmission-eigenvalue 
sequences with A^i and N2 levels (with Ni and A''2 con- 
strained to differ at most by one) is identical to a stag- 
gered superposition of two transmission-eigenvalue se- 
quences with A^i and N2 levels, which are correlated by 
the ordering requirement (jlSp . 



A. Constraints on the scattering matrix 

The presence of a lead-transposing symmetry immedi- 
ately results in the constraint r — r' , t — t' (when time- 
reversal symmetry is broken by a magnetic field, this can 
be achieved by an inversion symmetry but not by a re- 
flection symmetry). In order to further exploit the con- 
sequences of the symmetry, let us inspect a time-reversal 
symmetric system with a reflection symmetry, as shown 
in Fig. [TJe). As shown in Fig.[3ja,b), the desymmetrized 
versions are cut at the symmetry line, where they are 
equipped with Dirichlet or Neumann boundary condi- 
tions for wavefunctions of odd (— ) or even parity (-I-), 
respectively. Such wave functions are readily constructed 
starting from the original system when one chooses in- 
coming amplitudes of the form a'^^^ = zba*^-'^^ . The outgo- 
ing amplitudes are then given by b^^' = (r±t)a(^^ . Con- 
sequently, the scattering matrices of the desymmetrized 
systems are given by 



S±=r±t. 



(26) 



The desymmetrized systems only possess a single open- 
ing. In order to revert to the scattering matrix of the 
original system we invert Eq. (j26p . The transport in the 
original system is therefore described by the transmission 
matrix t = ^ (5+ — S- ) , which gives 



ttt = 1(2 



S-i^S'^ S — s I ) 



(27) 



The properties of this matrix — and especially, of its 
eigenvalues r„ — are not separable and depend on the in- 
terplay of both desymmetrized variants. 



B. Conventional random-matrix theory 



IV. LEAD-TRANSPOSING SYMMETRIES 

Systems with a lead-transposing symmetry require a 
separate treatment since the symmetry operation only 
commutes with the Hamiltonian, but not with the cur- 
rent operator (which changes its sign). In the presence 
of an applied bias, the symmetry operation exchanges 
the electronic source and drain reservoirs. An obvious 
symptom of this complication is the fact that the desym- 
metrized system only possesses a single lead (see Fig. 
[3]). Mathematically, the transmission matrix docs not 
assume a block structure but remains full. We will first 
adapt the concept of desymmetrization to derive the con- 
straints of the scattering matrix, and then turn to the 
joint probability density of the transmission eigenvalues 
in random-matrix theory. Just as in the previous section, 
we then focus on the orthogonal symmetry class {(3 = 1) 
and derive a closed expression for the joint probability 
density, which again assumes the form of a staggered level 
repulsion. 



Random-matrix ensembles for systems with lead- 
transposing symmetry can be obtained by assuming 
that the scattering matrices 5'+ and S- of the desym- 
metrized variants are statistically independent realiza- 
tions of the appropriate standard circular ensemble. The 
resulting ensembles are identical to those introduced by 
Baranger and Mello who based their considerations 
on a maximal-entropy principle. 

Earlier works have addressed isolated aspects of these 
ensembles, but not the complete transmission-eigenvalue 
statistics. For instance, it has been observed that a lead- 
transposing symmetry increases the conductance fluctu- 
ations |8| , |lll| but eliminates the weak-localization cor- 
rection For large N, the conductance fluctuations 
double, just as is the case for lead-preserving symmetries. 
We now provide a complete explanation of these observa- 
tions on the basis of the joint probability density of the 
transmission eigenvalues. 

The starting point of these considerations is the rela- 
tion 

T„ = sin2(e„/2) = i(l-cose„) (28) 
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between the transmission eigenvalues T„ and the eigen- 
phases 0„ of the unitary matrix Q = 5*+ 5*1, which fol- 
lows from Eq. ([27]). As illustrated in Fig. EJCjd), the 
statistics of transmission eigenvalues is hence directly im- 
posed by the statistics of the real parts cos8„ of the 
unimodular eigenvalues e'®" of Q. 

In random-matrix theory, the eigenphases 0„ follow 
the statistics of the associated circular ensemble. This 
is evident for the unitary ensemble {(3 — 2), which is in- 
variant under the multiplication of an arbitrary fixed ma- 
trix (it hence suffices, e.g., to assume that 5*+ is random 
while S*- is fixed, or vice versa). In the orthogonal case 
{(3 = 1), the unitary transformation Q' ~ S_^^^QS^^^ — 

SZ^^^S+SZ^^'^ results in a symmetric matrix with identi- 
cal eigenvalues. Their circular statistics then follows from 
the fact that the circular orthogonal ensemble is invariant 
under the symmetric involution with any fixed symmet- 
ric matrix (here, S_ ). The same transformation also 
succeeds in the case of self-dual matrices {(3 = 4). 

Because of the uniform distribution of eigenphases in 
the circular ensemble the one-point probability den- 
sity P{Tn) is given by Eq. Q fo r any finite N (i.e., not 
only in the limit N ^ oo) [llj. The joint probability 
density of the eigenphases 6„ is given by [1] 



^^e({e„}) « n 



cr,„, sm ■ 



e™ - e 



(29) 



Here we ordered the eigenphases by their moduli, 

< lOil < leal < leal < . . . < leArl < tt, (30) 

and denoted (t„ = sgne„. Since Eq. (I^S)) does not dis- 
criminate the sign of e„ we proceed to the distribution 
of the moduli 6'„ = |e„|. 



Pe{{On}) - ^e({fT„0„}). 



With the help of the relations 



sin(&„/2) = V?^, cos(0„/2) = - Tn, 
and also accounting for the Jacobian 
dOr,. 1 



dTn V7^„(i-r„)' 

this yields the joint probability density [ll| 
1 



(31) 



(32) 



(33) 



{cr„} m>n 



(34) 



This expression is symmetric under the replacement 
r„ — !■ 1 — T„, which explains the absence of weak- 
localization corrections to the conductance. Moreover, 



transmission eigenvalues do not repel each other when 
CT„ = —(Tm, i.e., when the underlying eigenphases e„ lie 
on the opposite (upper and lower) arcs of the unit cir- 
cle [see again Fig. Ellc)]. As the sets of eigenphases on 
both arcs is only weakly cross-correlated, this explains 
the doubling of the conductance fluctuations for large N. 



C. Staggered level repulsion for P — 1 

While the general conclusions of the previous section 
can be drawn for any (3, it should be noted that Eq. 
P4p still involves a combinatorial sum, and hence is sim- 
ilar in status as expression (jlip for systems with a lead- 
preserving symmetry. We now show that a much more 
detailed insight is possible for the orthogonal symmetry 
class {(3 — 1), where the combinatorial sum in Eq. ([M]) 
can be carried out explicitly (see below). The result- 
ing statistics again assume the form of a staggered level 
repulsion, but are not identical to Eq. (which was 
derived from the superposition of two independent level 
sequences): For N an odd integer, we find 



■m>n, both odd 



z odd VTiil-T^ 



(35a) 



m>n, both even 



while for even 

Pi{Tn}) (X 



n (T™-T„)n;4 



m>n, both odd 



/ odd 



n m--^") n -7^ 



m>n, both even 



/ even 



(35b) 



Similar to Eq. (|14p . the joint probability density again 
separates into two factors, each involving only every sec- 
ond eigenvalue. In particular, neighboring levels are not 
prohibited to approach each other closely, and statisti- 
cal fluctuations of observables are enhanced, as has been 
earlier observed for the conductance and the Fano fac- 
tor [1, [13, [ni . The correlation between the two level 
sequences is again imposed only indirectly by the require- 
ment that the sequences are staggered. This ordering re- 
quirement is independent of the parity of the wave func- 
tion - indeed, in the present case, parity is not well de- 
fined as the transmission eigenvalues arise from the com- 
bined properties of 5*4- and . 

In order to derive Eq. ([55]) . let us first inspect Eq. ([^ . 
Because of the ordering pO]) . each factor cr„i appears to — 
1 times, and therefore 



Pe{{On}) U^ili 

I even m>n 



sm ■ 



e„ 



(36) 



We next pass over to the joint distribution of moduli 
(I3ip . In order to evaluate the combinatorial sum over the 
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an we express the factor of sine functions in Eq. ((36)) as 
a Vandermonde determinant, 

Yl sin ®"'^ ^" ^ i-^ f(^-^y'detB{{a,,en}), (37) 

where -Bm;({0n}) = exp(z6mO! ™ = 1, 2, 3, . . . , iV, while 
the index / runs in integer steps from —{N — l)/2 to 
(iV — l)/2. The muhiUnearity of the determinant then 
yields 

Pe{{en}) oc (-z/2)^(^-i)/2 det C, (38) 

where Cmi = 2cos{9ml) for odd to and Cmi = 2^sin(^^,„Z) 
for even to. 

For every ^ > we now add the Zth column in C to the 
— Zth column, which cancels all sine terms in the latter 
columns. The determinant det C — det D det E then fac- 
torizes, where Dmi — coa 0m.l, rn odd, and Emi = sui0„il, 
TO even. If N is even, the index / is now restricted to 
I = 1/2, 3/2,..., (A^- l)/2. For odd N, this index is 
restricted to ? = 0, 1, 2, . . . , (A^ - l)/2 for the matrix D, 
and to / = 1, 2, . . . , (TV - l)/2 for the matrix E. 

For odd N we can write Dmi as a polynomial of de- 
gree / in cos 9m, and E^i as sin^^ times a polynomial of 
degree Z — 1 in cos 9m ■ We only need to keep the high- 
est monomial, as the other terms are linear combinations 
of the rows of lower index I. This leaves us again with 
Vandermonde determinants, 

detD (X (cos6'„ — cos^m), (39) 

m>n, both odd 

deti; cx Y\_ shi6'; Y\_ (cos6'„ - cos6',„). (40) 

I even m>n, both even 

For even N, the index I is half-integer, and the elements 
of D can now be written as cos(0m/2) times a polynomial 
in cos{9m), while those of E can be written as sin(0m/2) 
times such a polynomial. This yields 

detD cx ]J cos(e',/2) J]^ (cos 6'„ - cos 6'„j), (41) 

I odd m>n, both odd 

detE cx Yl sin(6'i/2) J| (cos 6l„ - cos 6l„0. (42) 

I even m>n, both even 

The joint probability density (|35p follows by transforming 
from 9n to T„, where the Jacobian is given by Eq. p3p . 
while the factors in the expressions for D and E can 
be rewritten with the help of Eq. (|32p and the relation 

C0S6I„ - C0s6'm = 2{Tm ~Tn). 

D. Large- A'^ asymptotics 

It is natural to ask whether the similarity of Eq. ([55]) 
to Eq. (|14p indicates a possible interpretation as a super- 
position of two independent level sequences [from which 
Eq. was derived]. In Eq. ([55]) . however, this inter- 
pretation is prevented by the different one-point weight 



terms associated to the even and odd indexed eigenval- 
ues. A symptom of this difference is the fact that Eq. 
([T4|) implies finite- iV weak-localization corrections to the 
conductance, while Eq. (|35p delivers the absence of such 
corrections, in agreement with the general conclusions in 
Sec. IIVBI Hence, the statistics of systems with a lead- 
transposing and a lead-preserving symmetry (with (3 = 1) 
only find a common ground when both are interpreted as 
a staggered level sequence. 

For the case of a lead-preserving symmetry, the frame- 
work of superpositions of independent level sequences 
of course provides a powerful tool for the derivation of 
low-point correlation functions and local statistics [such 
as the two-point correlation function, or the level spac- 
ing distribution (fT3|) ]. We now argue that in the limit 
iV — > oo, this framework can also be adopted for systems 
with a lead-transposing symmetry. 

In this limit, the transmission eigenvalues form a quasi- 
continuum, and the asymptotical statistics follow from 
the formal analogy to the statistics of coordinates of 
a dense set of parallel line charges in one dimension 
(the Coulomb gas), which exhibit a logarithmic repul- 
sion [H, Q. In leading order, the weight terms enter the 
analysis of the statistical fluctuations only via the one- 
point function P{T): For fixed index n, the transmis- 
sion eigenvalue r„ are confined to a small neighborhood 
around a nominal equilibrium position r„ , which is given 

by the implicit equation n - 1/2 = Jq^" P{T)dT. Sub- 
sequently, the weight terms can be approximated by a 
constant (with all the T„ fixed to T„), while the fluc- 
tuations are exclusively governed by the level-repulsion 
factors of the joint probability distribution. As the level- 
repulsion factors arc identical in Eqs. and (|35p one 
concludes that the local statistics in both ensembles be- 
come indistinguishable in the limit oi N oo. 

We therefore obtain the following remarkable result of 
purely statistical origin: For a lead-transposing symme- 
try, as N is sent to infinity the local statistics (embodied 
in low-point correlation functions) converges to that of a 
superposition of two independent level sequences. This 
is the case even though a classification of transmission 
eigenvalues by parity is not possible. In particular, we ar- 
rive at the prediction that in this limit, the level-spacing 
distribution is well approximated by Eq. p^ . 

V. NUMERICAL INVESTIGATIONS 

For the three standard Dyson ensembles of random- 
matrix theory, the joint probability density ([S]) manifests 
the celebrated repulsion between neighboring eigenval- 
ues, since the probability to find two closely spaced ad- 
jacent eigenvalues is suppressed as [Tn+i — Tn)^ . In con- 
trast, the joint densities (fH)) and ([55]) (both derived for 
(3=1) describe sequences of reduced stiffness, where only 
every second level is subject to mutual level repulsion. As 
argued before, as long as N takes on moderate values, 
the latter joint densities imply quantitative differences in 
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FIG. 4: (Color online) Probability density P(s) of 
transmission-eigenvalue spacings for systems with a lead- 
transposing symmetry, obtained from 10'' random matrices 
with TV = 50. Smooth curves: Spacing probability den- 
sity (|13|) for superpositions of eigenvalues of two independent 
sequences from the standard circular ensembles. The inset 
shows the Wigner distributions (|10p from standard random- 
matrix theory and the Poisson distribution (|9]) for uncorre- 
lated eigenvalues. 



the transmission eigenvalue statistics for lead-preserving 
and lead-transposing symmetries, while for large N these 
statistics should converge onto each other. 

In this section we illustrate the differences and simi- 
larities between these scenarios for all three main sym- 
metry classes {(3 = 1, 2, 4) via numerical sampling of the 
random-matrix ensembles, and also compare to realis- 
tic model systems. For convenient characterization of 
the eigenvalue repulsion we employ the nearest-neighbor 
spacing distribution P{s), as well as spacing distribu- 
tions to more distant neighbors. As we will see, the lo- 
cal statistics of systems with a lead-transposing symme- 
try actually show a much weaker N dependence than for 
systems with a lead- preserving symmetry. This feature 
could be anticipated by (but also goes beyond) the ab- 
sence of weak localization corrections in the one-point 
function (discussed in Sec. IIVB|) . 



N=4, |3=1 N=100, P=l 




N=4, |3=2 N=100, P=2 




N=4, |3=4 N=100, |3=4 




FIG. 5: (Color online) Probability densities of spacings s„ 
to the first, second and third neighboring transmission eigen- 
value for the random-matrix ensembles of systems with a lead- 
transposing symmetry (solid curves) or a lead-preserving sym- 
metry (dashed curves). In the left panels the number of trans- 
port channels A'' = 4, while in the right panels N = 100. Top 
panels: orthogonal symmetry class (/3 = 1). Middle panels: 
unitary symmetry class (/3 = 2). Bottom panels: symplec- 
tic symmetry class (/3 = 4). For each ensemble, the results 
represent a sample of 10'* realizations. 



A. Random-matrix theory 

We start with the characterization of the level statistics 
within the various random-matrix ensembles. Let us first 
consider the case of a lead-transposing symmetry with a 
relatively large number of transport channels, for which 
we expect that the local statistics is close to that of a su- 
perposition of two independent level sequences. Starting 
point of the numerical computations is Eq. I|27p . where 
the matrix Q = S+Sl_ is drawn from the appropriate 
Dyson ensemble. In order to obtain the nearest-neighbor 
spacing distribution P(s), we unfold the eigenvalue se- 
quences to a mean local spacing s = 1 [1, Q . Figure [3] 
shows the resulting spacing distributions for A'^ = 50. For 



this large number of open channels we find that the nu- 
merical histograms indeed match the predictions from the 
superposition of two independent level sequences [solid 
curves; see Eq. ([T^ ]. 

For comparison, the inset in Fig. 3] shows the stan- 
dard Wigner distributions (fTU)) . as well as the Poisson 
distribution ([9]). In the Poisson distribution the eigen- 
value spacing density is maximal at s = 0; for larger s 
the probability density decreases monotonically. For the 
Wigner distributions the most likely eigenvalue spacing 
occurs at a finite value of s; for s — > 0, the distributions 
decay algebraically cx s^, while for s — s- oo they decay as a 
Gaussian. The distributions in the main panel combine 
the partial absence of level repulsion for small s [with 
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quantum billiards open kicked rotators 
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FIG. 6: (Color online) (a) Nearest-neighbor spacing distri- 
bution P{s) for the lead-asymmetric stadium billiard of Fig. 
[TJb), averaged over energies in the range A'^ = 5 — 14, and 
the lead-transposing symmetric stadium billiard of Fig. Wie) , 
with A'' = 5, 6. (b) The same for open quantum kicked ro- 
tators with TV — 12. In both panels, the solid curves show 
the Wigner distribution (|10p with (3 = 1 and the prediction 
of random-matrix theory for systems with a lead-transposing 
symmetry [which can be safely approximated by Eq. (|13a|) ]. 



P(s = 0) = 1/2] with the Gaussian decay of the Wigner 
distributions for large s. 

For large N , virtually identical results are obtained 
for the conventional case of a lead-preserving symmetry. 
This is demonstrated in detail in Fig.O which also shows 
the spacing distributions to the second and third-nearest 
neighbor. Here, solid curves are for a lead-transposing 
symmetry, and dashed curves are for a lead-preserving 
symmetry (corresponding to a superposition of indepen- 
dent level sequences from the appropriate Dyson ensem- 
ble). For N = 100 (right panels), dashed and solid curves 
lie on top of each other and are practically indistinguish- 
able. This clearly supports the convergence of the local 
statistics of both cases for large N . 

The left panels in Fig. [5] show the level-spacing dis- 
tributions for = 4. In this case, the results for 
a lead-transposing symmetry are distinctively different 
from those for a lead-preserving symmetry. Interest- 
ingly, the nearest-neighbor spacing distribution for a 
lead-transposing symmetry is very similar for small and 
large N\ the distribution for = 4 is already well ap- 
proximated by Eq. ([13]). In comparison, the nearest- 
neighbor spacing distribution for a lead-preserving sym- 
metry shows a much stronger A^-dependence. 



B. Comparison to model systems 

In order to validate that realistic quantum systems can 
indeed be described by random matrix theory (on which 
all previous considerations are based), we compare our 
predictions with numerical results for such systems. In 
particular, we present results of numerical computations 



for quantum billiards, which model a lateral quantum 
dot, and for the open kicked rotator, which is based on 
an efficient quantum map. We focus on systems in the 
orthogonal symmetry class (/3 = 1) and contrast sys- 
tems with a lead-transposing symmetry to systems with- 
out any spatial symmetry. 

The quantum billiards are derived from the stadium 
geometry, with leads positioned to either break or con- 
serve the reflection symmetry about the vertical center 
line [see Figs. [Ub,e)]. The computations are performed 
using a modular recursive Green's function method [H, 
[Tsj . with energies that permit 5 < < 14 open chan- 
nels in each of the two leads. As shown in Fig. [6ja), the 
eigenvalue spacing distribution agrees well with the pre- 
dictions of random-matrix theory, both in presence and 
in absence of a lead-transposing symmetry. 

The open quantum kicked rotator [3, [l^, [l^ is defined 
by the scattering matrix 

S = P[e-'^ - F{1 - P^P)]-^FP^, (43) 

where e is the quasi-energy, 

Fn„^ = (iM)-i/2eif Tr+™« t) (44) 

is the M X M-dimensional Floquet operator of the kicked 
rotator, and P is an 2N x Ajf-dimensional matrix which 
projects the internal Hilbert space onto the openings. 
We assume that M is even and M N . The reflec- 
tion symmetry of the closed system is manifested in the 
symmetry F„m = FM^n,M~m, and the lead-transposing 
symmetry of the open system is present when in addition 

Pnm = P2N-n,M-m- 

Figure [Hlb) shows the spacing distributions obtained 
for kicked rotators with symmetrical and asymmetrical 
lead placement and A^ = 12. The data represents 6600 re- 
alizations which are generated by varying the quasienergy 
e e [0, 27r), the kicking strength K G [10, 15], and the in- 
ternal dimension M e [498,502]. Again, we find good 
agreement with random-matrix theory, including the re- 
duced eigenvalue repulsion in the lead-symmetric case. 

The results in this section reveal clear signatures of 
staggered level repulsion in realistic systems with a lead- 
transposing symmetry (and P = 1). It is worth emphasiz- 
ing that the applicability of this statistical concept [em- 
bodied in the random-matrix results Eq. ([35])] does not 
rely on any pre- or postprocessing or -selection of the 
transmission eigenvalues in the model systems (as there 
is no intrinsic property of the transmission eigenvalues 
or their associated scattering wave functions - such as a 
parity - that could be used to divide these eigenvalues 
into two sets). 

VI. SUMMARY AND CONCLUSIONS 

We analyzed the transport in open systems with a lead- 
transposing or a lead-preserving symmetry via the com- 
plete joint probability density of transmission eigenval- 
ues, obtained in random-matrix theory. 
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For a Icad-prcscrving symmetry, the standard concept 
of desymmetrization reduces the problem to the inves- 
tigation of independent non-symmetric variants of the 
system. For a lead-transposing symmetry, however, the 
transport characteristics only arise as a collective prop- 
erty of the symmetry-reduced variants of the system. We 
still found that both types of symmetry result in a similar 
reduction of level repulsion, so that transmission eigen- 
values can approach each other closely. For a large num- 
ber of transport channels N, the local eigenvalue statis- 
tics for both types of symmetry indeed become indistin- 
guishable. 

Our main analytical results concern a detailed expla- 
nation of these features for systems which also exhibit 
time-reversal and spin-rotation invariance (the orthogo- 
nal symmetry class, with symmetry index /? = 1). In 
this case, the transmission eigenvalue statistics of sys- 
tems with a lead-transposing or lead-preserving symme- 
try find a common natural interpretation as a staggered 
superposition of two independent level sequences. In such 
a superposition the eigenvalues alternate between the se- 
quences when they are ordered by magnitude. The joint 
probability densities for the two types of symmetry only 
differ in one-point weight factors. For lead-preserving 



symmetries these weight factor incorporate l/N correc- 
tions for quantities such as the ensemble-averaged con- 
ductance, while these corrections are absent for a lead- 
transposing symmetry. This results in differences of the 
local eigenvalue statistics when N is small, but becomes 
insignificant when A'" is large. 

While we concentrated on systems with discrete spatial 
symmetries, our results can also be applied for discrete 
symmetries of different origin (e.g., arising from internal 
degrees of freedom) that yield equivalent constraints on 
the scattering matrix. 
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